Fitting the Linear Model of Coregionalization by Generalized Least Squares   总被引：2，自引：0，他引：2  

Bernard Pelletier  Pierre Dutilleul  Guillaume Larocque  James W. Fyles 《Mathematical Geology》2004,36(3):323-343

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.  相似文献   

Conditions of pocket formation in the Oktyabrskaya tourmaline-rich gem pegmatite (the Malkhan field, Central Transbaikalia, Russia)     

Igor S. Peretyazhko  Victor Ye. Zagorsky  Sergey Z. Smirnov  Mikhail Y. Mikhailov 《Chemical Geology》2004,210(1-4):91-111

Coexisting melt (MI), fluid-melt (FMI) and fluid (FI) inclusions in quartz from the Oktaybrskaya pegmatite, central Transbaikalia, have been studied and the thermodynamic modeling of PVTX-properties of aqueous orthoboric-acid fluids has been carried out to define the conditions of pocket formation. At room temperature, FMI in early pocket quartz and in quartz from the coarse-grained quartz–oligoclase host pegmatite contain crystalline aggregates and an orthoboric-acid fluid. The portion of FMI in inclusion assemblages decreases and the volume of fluid in inclusions increases from the early to the late growth zones in the pocket quartz. No FMI have been found in the late growth zones. Significant variations of solid/fluid ratios in the neighboring FMI result from heterogeneous entrapment of coexisting melts and fluids by a host mineral. Raman spectroscopy, SEM EDS and EMPA indicate that the crystalline aggregates in FMI are dominated by mica minerals of the boron-rich muscovite–nanpingite CsAl₂[AlSi₃O₁₀](OH,F)₂ series as well as lepidolite. Topaz, quartz, potassium feldspar and several unidentified minerals occur in much lower amounts. Fluid isolations in FMI and FI have similar total salinity (4–8 wt.% NaCl eq.) and H₃BO₃ contents (12–16 wt.%). The melt inclusions in host-pegmatite quartz homogenize at 570–600 °C. The silicate crystalline aggregates in large inclusions in pocket quartz completely melt at 615 °C. However, even after those inclusions were significantly overheated at 650±10 °C and 2.5 kbar during 24 h they remained non-homogeneous and displayed two types: (i) glass+unmelted crystals and (ii) fluid+glass. The FMI glasses contain 1.94–2.73 wt.% F, 2.51 wt.% B₂O₃, 3.64–5.20 wt.% Cs₂O, 0.54 wt.% Li₂O, 0.57 wt.% Ta₂O₅, 0.10 wt.% Nb₂O₅, 0.12 wt.% BeO. The H₂O content of the glass could exceed 12 wt.%. Such compositions suggest that the residual melts of the latest magmatic stage were strongly enriched in H₂O, B, F, Cs and contained elevated concentrations of Li, Be, Ta, and Nb. FMI microthermometry showed that those melts could have crystallized at 615–550 °C.

Crystallization of quartz–feldspar pegmatite matrix leads to the formation of H₂O-, B- and F-enriched residual melts and associated fluids (prototypes of pockets). Fluids of different compositions and residual melts of different liquidus–solidus P–T-conditions would form pockets with various internal fluid pressures. During crystallization, those melts release more aqueous fluids resulting in a further increase of the fluid pressure in pockets. A significant overpressure and a possible pressure gradient between the neighboring pockets would induce fracturing of pockets and “fluid explosions”. The fracturing commonly results in the crushing of pocket walls, formation of new fractures connecting adjacent pockets, heterogenization and mixing of pocket fluids. Such newly formed fluids would interact with a primary pegmatite matrix along the fractures and cause autometasomatic alteration, recrystallization, leaching and formation of “primary–secondary” pockets.  相似文献   

天祝-古浪5.4级地震前震源区应力场的短临变化特征   总被引：1，自引：0，他引：1  

荣代潞  张元生  李亚荣 《高原地震》2000,(2)

利用矩张量反演方法反演了 1996年 6月 1日天祝 -古浪 5 .4级地震序列的震源机制 ,并与该区小震综合断层面解进行了比较。在此基础上讨论了该地震前后震源区应力场的变化特征。  相似文献   

对矿物包裹体溶液性质的新认识——兼论含矿溶液是酸性溶液的地质证据     

陈福  龙洪波  董丽敏 《矿物学报》2000,(1)

对金属矿床成矿溶液的性质 ,特别是 pH值 (酸碱度 )的确定 ,是认识含矿溶液的成因、成矿条件和成矿机理的关键。目前对成矿溶液性质的研究主要局限在对围岩蚀变矿物包裹体的测定方面 ,并得出近中性热卤水成矿的结论。作者认为 ,多数围岩蚀变矿物的包裹体溶液不能代表原成矿溶液的性质 ,“近中性的热卤水成矿”的认识有误 ,应是酸性水形成的含矿溶液成矿 ,并从矿物包裹体溶液的特征、Roedder的高金属含量的矿物包裹体溶液的发现、氢氧同位素资料、矿物包裹体中高卤水溶液的起源以及成矿溶液的演化机理等方面论述了酸性金属含矿溶液的成因及成矿演化过程。  相似文献   

Global Bifurcations of Periodic solutions of the Hill Lunar Problem*     

A.J. Maciejewski  S.M. Rybicki 《Celestial Mechanics and Dynamical Astronomy》2001,81(4):279-297

  相似文献   

Physics of fracturing and seismic energy release: A review     

S. J. Gibowicz 《Pure and Applied Geophysics》1986,124(4-5):611-658

This review was prepared as an opening lecture for the International Symposium on Physics of Fracturing and Seismic Energy Release, held at the Castle of Liblice near Prague from October 28 to November 1, 1985, and organized by the Geophysical Institute of the Czechoslovak Academy of Sciences. The review attempts to classify and synthesize results of recent studies in fracture mechanics and earthquake source physics. The following topics are discussed: recent developments in fracture mechanics, earthquake source modeling, heterogeneous fault planes, foreshocks and aftershocks, faults and fractals, and moment tensor solutions. This rather broad range of topics prevents presentation of a complete list of all relevant works, though over one hundred and seventy references are cited.  相似文献   

Global Bifurcations of Periodic Solutions of the Restricted Three Body Problem     

A. J. Maciejewski  S. M. Rybicki 《Celestial Mechanics and Dynamical Astronomy》2004,88(3):293-324

We describe global bifurcations from the libration points of non-stationary periodic solutions of the restricted three body problem. We show that the only admissible continua of non-stationary periodic solutions of the planar restricted three body problem, bifurcating from the libration points, can be the short-period families bifurcating from the Lagrange equilibria L ₄, L ₅. We classify admissible continua and show that there are possible exactly six admissible continua of non-stationary periodic solutions of the planar restricted three body problem. We also characterize admissible continua of non-stationary periodic solutions of the spatial restricted three body problem. Moreover, we combine our results with the Déprit and Henrard conjectures (see [8]), concerning families of periodic solutions of the planar restricted three body problem, and show that they can be formulated in a stronger way. As the main tool we use degree theory for SO(2)-equivariant gradient maps defined by the second author in [25].This revised version was published online in October 2005 with corrections to the Cover Date.  相似文献   

The general solution of the HENON–HEILES problem     

K.?E.?PapadakisEmail author  C.?L.?Goudas  G.?A.?Katsiaris 《Astrophysics and Space Science》2005,295(3):375-396

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.  相似文献   

The evolution characteristics of positive and negative earthquakes before strong quakes in Kanto area and the Hyugo earthquake in Japan     

王裕仓  宋治平  细野耕司  陈学中  尹祥础 《地震学报(英文版)》1998,11(1):35-42

ＴｈｅｅｖｏｌｕｔｉｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｐｏｓｉｔｉｖｅａｎｄｎｅｇａｔｉｖｅｅａｒｔｈｑｕａｋｅｓｂｅｆｏｒｅｓｔｒｏｎｇｑｕａｋｅｓｉｎＫａｎｔｏａｒｅａａｎｄｔｈｅＨｙｕｇｏｅａｒｔｈｑｕａｋｅｉｎＪａｐａｎＹＵＣＡＮＧＷ...  相似文献   

端承桩在层状地基波动与辐射阻尼影响下的横向地震反应分析   总被引：1，自引：0，他引：1  

刘宗贤  李玉亭 《东北地震研究》1998,14(2):43-52

通过对层状地层中端承桩的分析,建立了合理的几何,数学及力学模型。考虑层状地基ＳＨ与ＳＶ波波动及非线性辐射阻尼影响,应用Ｄｉｒａｃ广义函数及Ｗｉｎｋｌｅｒ地基理论计入各种动力,抗力,地震力及内力等,建立地基波动方程,桩振动方程,地基与桩耦合振动方程。给出了柯西型函数的自由振动解。利用非经典理论正交条件和卷积定理给出了强迫振动地震反应解析解,算例结果充分表明了动力特性及振动规律非常正确,充分体现了本文  相似文献